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Abstract 

This paper presents a recursive algorithm for the deadbeat predictive control that brings the 
output response to rest after a few time steps. The main idea is to put together the system 
identification process and the deadbeat control design into a unified step. It starts with 
reformulating the conventional multi-step output-prediction equation to explicitly include 
the coefficient matrices to weight past input and output time histories for computation of 
feedback control force. The formulation thus derived satisfies simultaneously the system 
identification and the deadbeat control requirements. As soon as the coefficient matrices arc 
identified satisfying the output prediction equation, no further work is required to design 
a deadbeat controller. The method can be implemented recursively just as any typical 

recursive system identification techniques. 
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1 Introduction 


There arc many interesting technical problems in the area of controlled aerospace structures 
that NASA researchers are trying to solve. These problems, for example, include acous- 
tic noise reduction, flow control, ride quality control, flexible spacecraft attitude control 
and vibration suppression. Active or passive control for a dynamic system is not a new 
subject. Many control techniques 1-5 are available today and ready to be used for applica- 
tion to these interesting problems. Some of the techniques arc the quadratic optimization 
technique 1-2 , the pole placement technique 3 , the virtual passive technique 4 , the energy dis- 
sipation technique 5 , and the adaptive control technique . 6-18 Some researchers prefer to work 
in the frequency domain using the frequency response functions (FRF) while others use the 
state-space model (SSM) in the time-domain to design controllers. The model-based tech- 
niques need a mathematical model (FRF of SSM) within a certain level of accuracy to design 
a controller. Except for a few simple cases, system identification must be involved in the 
design process to verify the open-loop model and the closed-loop design as well. As a result, 
it may take considerable time to iterate the design process until performance requirements 
arc met. For the systems with minimum uncertainties, the iteration procedure would not 
bother the control engineers, as long as a satisfactory control design can be found. 

For systems with unknown disturbances and considerable uncertainties, the controller 
must be able to adapt the unknown changes in real time. Adaptive control techniques 
are developed for this purpose. The approach is to adjust the control gains to reflect the 
system changes so as to continuously check and meet the performance requirements. Most 
adaptive control techniques require the controlled system to be minimum phase 19-25 in the 
sense that all the system transmission zeros are stable. The minimum-phase system in 
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the continuous-time domain does not guarantee its minimum phase in the discrete-time 
domain. In practice, only a few structural systems in the discrete-time domain arc minimum 
phase. Predictive controller designs 11-18 were developed to particularly address the non- 
minimum-phase problems with the hope that they can be implemented in real time. Two 
indirect techniques and one direct technique were derived in Ref. 18 using the concept of 
deadbeat predictive control law that brought the output response to rest after a few finite 
time steps. The indirect techniques require identification of coefficient matrices of a finite- 
difference model representing the controlled system. The deadbeat predictive controllers arc 
then computed using the identified coefficient matrices. Note that the identified matrices 
minimize the output error between the estimated and real outputs. The direct technique 
computes the deadbeat predictive controller directly from input and output data without 
explicitly identifying the system parameters. However, it requires to minimize the output 
error first and then perform a Hankel-like 18 matrix to calculate the control gains for the past 
input and past output signal. Since it takes time to invert a matrix, both direct and indirect 
algorithms have a drawback for real time application. Nevertheless, the direct algorithm did 
provide the fundamental framework for further development of a recursive technique for real 
time implementation. 

A new recursive technique is presented in this paper for the design of a deadbeat 
predictive controller. It uses the approach derived in Ref. 18 for the direct algorithm. The 
technique computes the gain matrices recursively and directly from input and output data in 
every sampling period. In addition, the recursive formula satisfies both system identification 
and deadbeat predictive control equations simultaneously. As a result, the design process is 
completed in such a way that there is no time delay between the identification step and the 
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control gain computation step. 

This paper begins with a brief introduction of multi-step output prediction. 18 The 
basic formulation for a deadbeat controller design is then derived giving the mathematical 
foundation of the recursive method. A recursive formula with computational steps is also 
included for real-time implementation. With a slight modification, the formula is extended 
to compute the feedforward gain for a measurable or predictable disturbance input. Finally, 
several numerical examples are given for illustration of the method. 

2 Multi-Step Output Prediction 18 

The input/ output relationship of a linear system is commonly described by a finite difference 
model. Given a system with r inputs and m outputs, the finite difference equation for the 
r x 1 input u(k ) and the m x 1 output y(k) at time k is 

y(k + 1) = a x y(k) + a 2 y(k - 1) H F a p y(k - p + 1) (1) 

+P 0 u(k) + Piu(k) + (3 2 u(k — 1) -I + f3 p u{k — p + 1) 

It simply means that the current output can be predicted by the past input and output 
time histories. The finite difference model is also often referred to as the ARX model where 
AR refers to the AutoRegressivc part and X refers to the eXogeneous part. The coefficient 
matrices, a* (* = 1, 2, • • • , p) of m x m and j3{ (i = 1,2,---, p) ofmxr are commonly referred 
to as the observer Markov parameters (OMP) or ARX parameters. The matrix /3 0 is the 
direct transmission term. 

Equation (1) produces the following multi-step output prediction 18 

y p (k + q) = T'u p+q (k) + B'u p (k -p) + A'y p (k - p) (2) 
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The quantity y p (k + q) represents the output vector with a total of p data points for each 
sensor from the time step k + q to k + q + p — 1, whereas y p (k — p) includes the p data 
from k — p to k — 1. Similarly, u q+p (k) has q 4- p input data points starting from the time 
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step k and u p (k — p) has p input data points from k — p. The matrix T' is formed from the 
parameters, /?o \ • i and (3q +p 1 (the pulse response sequence). 

The vector y p {k+q) in Eq. (2) consists of three terms. The first term is the input vector 
u q+p (k ) including inputs from time step k to k + q + p— 1. Relative to the same time k , the 
second and third terms, u p (k — p ) and y p (k — p ), are input and output vectors from time step 
k — p to k — 1, respectively. The future input vector w s (fc) is to be determined for feedback 
control. The matrices B' and A! may be computed from observer Markov parameters «, 
(i = 1, 2, • • • ,p) and /?* (i = 1, 2, • • • ,p), or directly from input and output data. 

3 Deadbeat Predictive Control Designs 

Several deadbeat control algorithms have been developed 18 using the multi-step output pre- 
diction, Eq.(2). Among these algorithms, the direct algorithm uses the input and output 
data directly without using a t ( i = 1, 2, • • • ,p) and f3 t (i = 1, 2, - • • ,p) to first compute B' 
and A! and then design a deadbeat predictive controller. The goal was to make the direct 
algorithm suitable for real-time implementation in the sense that the deadbeat controller 
may be updated at every sampling interval. Unfortunately, it involves a matrix inverse that 
is difficult, if not impossible, to compute it recursively. To overcome the computational 
difficulty, an alternative algorithm is developed in this section. 

Let T' be partitioned into two parts such that Eq.(2) becomes 

y p {k + q) = T 0 u p (k + q) + T c u q (k ) + B'u p (k - p) + A'y p (k - p) (5) 


6 



where 
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Both T 0 of pm x pr and % of pm x qr are formed from system pulse response (system Markov 
parameters). Note that, m is the number of outputs, p is the order of the ARX model, r is 
the number of inputs, and q is an integer. Given any input and output sequence u(k) and 
y(k), Eq. (5) must be satisfied and can be used for identifying coefficient matrices T 0 , %. A ' , 
and B' . 

Note that the matrix T c of pm x qr is a Hankcl-like matrix 18 which has rank n where 
n is the order of the system. For the ease where qr > pm > n, there are qr elements in u q (k) 
of qr x 1 with only n independent equations in Eq. (5). As a result, Eq. (5) provide multiple 
solutions for u q (k ) with the minimum-norm solution expressed by 


u q {k ) = T}y p {k + q) - TjT 0 u p (k + q) - T^B'u p (k - p) - T? A!y p (k - p) 
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or in matrix form 


u q (k) = [ -T C U' -T]B' T c t -TJT 0 


( 9 ) 


Vv( k - P ) 
u p (k - p) 

y P i k + q) 

U p (k + q) 

For the case where qr = pm,, Eq. (9) is unique. To simplify Eq. (9), define the following 
notations 


F r = 
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and 
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( 12 ) 


and v p (k + q) = 

where both F c and F 0 are qr x p(m + r ) matrices, and both v p (k - p) and v p (k + q) are 

(pm + pr) x 1 column vectors. Equation (9) thus becomes 

v p (k - p) ' 

Vp(k + q) 

Equation (12) is another form of finite difference model for system identification. For any 
given input and output data, there exists a set of F c and F 0 satisfying Eq. (12). Using Eq. 
(12) to develop a deadbeat controller is shown in the following. 

Let us assume that the input vector u q (k) is chosen such that 


u q (k) — F c Vp(k - p) 


( 13 ) 


F oVp (k + q)= \ T c t -TJT 0 


(14) 


In order to satisfy Eq. (12), the following equation must hold 

y P {k + q) 

u p (k + q) 

If Tj of qr x pm is full rank pm with qr > pm , and u p (k + q) is set to zero, then y p (k + q) 
becomes zero. As a result, the control action u q (k) computed from Eq. (13) is a deadbeat 
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controller which makes y p {k + q) to zero after q time steps. The integer q is thus called 
the deadbeat control horizon, whereas the integer p is commonly referred to as the system 
identification (or observer) horizon. The condition, qr > pm, means that the control action 
should not be faster than the state observation. That makes the physical sense. 

In view of Eq. (6), the first r rows of u q (k) is the input vector u(k) at time k. Define 
F c i and F 0 \ as the first r rows of F c and F a , respectively. The control action at u(k) should 
be 


u(k) = [the first r rows of F c ] v p (k — p) 
= F cA v p (k - p) 

where F c \ is the control gain matrix to be determined. 

The first r rows of Eq. (12) produces 

u(k) = (the first r rows 
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(16) 


v p (k - p) 
v p (k + q) 

Equation (16) indicates that the input u(k) is related to the past input sequence u(k — p) 
to u(k — 1) and output sequence y(k — p) to y(k — 1), and future input sequence u{k + q) 
to u(k + q + p — 1) and output sequence y(k + q) to y(k + q + p — 1). There is a total of q 
time steps gap from A; to k + q — 1. A different integer q produces a different set of coefficient 
matrices F c i and Foi that satisfies Eq. (16). 
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To solve Eq. (16), let us first form the following matrices 


U(k ) = u(k) u(k + 1) ••• u(N — p — q + 1) 

y(k-p) y(k-p+ 1) ••• y(N-2p-q+l) 
u(k - p) u(k-p + 1) ••• u(N - 2p - q + 1) 
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u(k + q + p — 1) u(k + q + p) ■■■ u(N ) 

where N is the data length used for estimation of coefficient matrices F cl and F ol . Application 
of Eq. (16) thus yields 

V p (k - p) 

V p (k + q) 

Let the integer N be chosen large enough such that the matrix U(k) of r x ( N — p — q — k + 2) 
has rank r, V p {k-p ) of p(r+ra) x (N-p-q-k+2) andV p (k+q) oip(r+m)x(N-p-q-k+2) 
have rank pr + n where n is the order of the system. Again, r means the number of inputs 
and m represents the number of outputs. Equation (18) produces the following least-squares 
solution 
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4 Recursive Least-Squares Algorithm 


There arc many recursive algorithms 2 available to solve the least-squares problem. The 

classical least-squares method is the most straightforward approach and is also the basis for 

the others. The classical recursive method is briefly described here. 

Equation (16) can be written in a compact matrix form 

v p (k - p) ' 
v p (k + q ) 

= Fv p (k - 1) (20) 


u(k) = F c i F a i 
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cl 
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u(k + 1) = F(k)v p (k) 


( 22 ) 

(23) 


Next, compute the following quantities 


Pp(k) 

= P p (k-l)[I-Vp(k)G p (k)} 

(24) 

F(k + 1) 

= F(k) + [u(k + 1) — u{k + 1)]G P (A:) 

(25) 


Equations (22) (25) constitute the fundamental Recursive Least-Squares (RLS) formulation 
for identifying the gain matrix F including F c \ for the deadbeat controller design and F ol for 
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the need of system identification. The initial values of P p (0) and Y p (l) can be either assigned 
or obtained by performing a small batch least-squares after gathering a sufficient number 
of data. If initial values are to be chosen, P p ( 0) and F(l) can be assigned as dl 2p{r+m) 
and 0 rX 2p( r +m)> respectively, where d is a large positive number. The positive constant d 
is the only parameter required for the initialization. The proper choice of d is based on 
practical experience. The initialization introduces a bias into the parameter estimate F'(k) 
produced by the recursive least-squares method. For large data lengths, the exact value of 
the initialization constant is not important. 

4.1 Computational Steps after initialization 

The computational steps for the recursive deadbeat control method are summarized in the 
following. 

1. Form the vector v p (k) as shown in Eq. (21) with the new input v(k + q +p - 1) as the 
last r -f m rows. 

2. Compute the gain vector G p (k) by inserting P p (k - 1) and v p (k) in Eq. (22). In this 
step, one should compute v^{k)P p {k - 1) first and then use the result to calculate 
[v T p {k)P p {k-\)]v p (k). 

3. Compute the estimated output u(k + 1) by substituting F(k) and v p (k ) into Eq. (23). 

4. Update P p (k — 1) to obtain P p (k) with v p (k) formed from the first step and G p (k ) 
computed from the second step. 

5. Update F(k ) to obtain F(k + 1) from Eq. (25) with the input signal u(k + 1), the 
estimated input u(k + 1), and the computed gain G p (k). 
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No matrix inverse is involved in these computational steps. Updating P p (k) and G p (k) 
takes more time than computing other quantities. The recursive procedure derived for up- 
dating the least-squares solution F(k) is very general in the sense that it is valid for any 
linear equation such as Eq. (19). 

5 Feedback and Feedforward for Disturbance Input 

In addition to the control input, there may be other disturbance inputs applied to the system. 
Some type of disturbances comes from the known sources that can be measured. Another 
type of disturbances is not known but its correlation is known. This section addresses 
the predictive feedback designs including feedforward from the disturbance inputs that are 
measurable or predictable.. 

With the disturbance input involved, the multi-step output prediction equation be- 
comes 

y p (k+g) = T c u cq (k)+Tju d{q+p) (k)+T 0 u cp (k+q)+B' c u cp (k-p)+B d u dp (k-p)+A'y p {k-p) (26) 
where 
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The subscripts c and d signify the quantities related to the control force and the disturbance 
force, respectively. The form of the matrix T d associated with the disturbances u d is similar 
to T' defined in Eq. 4. The matrix T7 is a pm x qr c matrix where r c is the number of control 
inputs and T d is a pm x qr d matrix where r d is the number of disturbance inputs. The 
forms of B' c and J3' d are also similar but corresponding to different type of forces. Note that 
T C ,T 0 , and B' c are quantities associated with the control force u c (k). 

A similar equation to Eq. (8) can thus be derived as 

u cq (k) = T}y v {k + q) - T c ] T oUcp {k + q) - TjTju d{(]+p) (k) 

—Tj A' y p (k -p)- TjB'cUcpik - p) - TjB' d u dp (k - p) (28) 

or in matrix form 

’ y P {k - p) 

Ucp{k - p) 

u cq (k)=[-TjA' -TjB' c -TjB d T c t -TjT a -7?T' ] Ucd{k ~ p) ( 29 ) 

Vp{k + q) 

u C p(k + 5) 

. u d(q+p){k ) 

Define the following notations 

F' c = [ -TJA' -T<B' e -pB' d ] , Fi = [ Tj -T}T a ] , and F' d = -pp (30) 

and 

Vvi k ~ P ) 

u cp (k - p) 

Ucd(k - p) 

where F' c is a qr c x p(m 4- r c + r d ) matrix but F' co is a qr c x j>(rn I r c ) matrix. The quantity 
v p (k -p) is a p(m + r c + r d ) x 1 column vector whereas v^k -f q) is a p(m + r c ) x 1 column 


and Vcp(k + q) = 


y P ( k + q) 

Ucp(k + q) 


v p (k - p) = 
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vector. Equation (29) thus becomes 

v p (k - p) 

Ucq(k) = [ F' c F' co F' d ) Vcp{k + q) (32) 

u d(q+p ){ k ) 

For any given input and output data, there exists a set of F' c , F' co) and F' d satisfying Eq. (32). 
Let us assume that the input vector u C(j {k) is chosen such that 

Ucq{k ) = F' c v p (k - p) (33) 

From Eq. (32), the output after q time steps is then governed vy 

Ko v cp{k + q) + F' d u d[q+p) {k ) = 0 

or, from Eq. (29), 

Vv P ( k + q) = T?T 0 u cxl {k + q) + T^T d u d{q+p) (k) (34) 

Here we have assumed that the disturbance is measurable. If the disturbance is not pre- 
dictable or measurable, Eq. (33) is not valid. Equation (34) indicates that the output vector 
from time k + q is generated by the control vector from k + q and the disturbance vector from 
k. If the disturbance Ud{k) such as the random signal is not predictable, then we cannot use 
the disturbance signal for a feedforward design for the control force u(k) at time k. As a 
result, the feedforward design included in the control law, Eq. (33), is the only way that can 
be implemented, because it uses only the disturbance signal before time k. From Eq. (33), 
the control action u(k ) at time k is 

u(k) — (the first r rows of F' c ) v p (k — p ) 

= F' rA v p (k — p) (35) 
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where F' cl is the control gain matrices to be determined. 
The first r rows of Eq. (32) produces 
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A similar equation to Eq. (19) can then be obtained 
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Equation (37) can also be computed recursively. 


6 Numerical Example 


A simple spring-mass-damper system is used to illustrate various controllers. Several different 
cases will be discussed ranging from single-input/singlc-output to multi-input/multi-output. 
First, the noise-free case is shown and then the case with additive measurement noise is 
discussed. 

Consider a three-degree-of-freedom spring-mass-damper system 


Mid + =.w + Kw = u 


where 
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and rrii, ki, Q, i = 1, 2, 3 are the mass, spring stiffness, and damping coefficients, respectively. 
For this system, the order of the equivalent state-state representation is 6 {n = 6). The 
control force applied to each mass is denoted by Ui,i = 1,2,3. The variables w , , i — 1,2,3 
are the positions of the three masses measured from their equilibrium positions. In the 
simulation, mi = m2 = m 3 = IKg, ki = k 2 = k 3 = 1,000 N/m, (i = C2 = C3 = 0.1AT — 
sec/m. The system is sampled at 50 Hz (At = 0.02sec.). Let the measurements fji be the 
accelerations of the three masses, y t = id j, i = 1, 2, 3 . 

Let us consider a singlc-control-input and single-output case where the control input 
to the system is the force on the first mass (i.c., u c = U\), and the output is the acceleration 
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of the third mass (i.e., y — u> 3 ) (non-collocated actuator-sensor). Therefore, the smallest 
order of the ARX model p is 6 corresponding to a deadbeat observer, and the smallest value 
for q is also 6 corresponding to a deadbeat controller which will bring the entire system to 
rest in exactly 6 time steps. Note that this is a non-minimum phase system. 

Consider the case where the controller is computed with q = 6. Let the initial guess 
for P p (0) and F(l) shown in Eqs. (22 ) and (23) be 1000/ 2 4 X 24 and 0i x2 4, respectively. The 
input signal is a sequence of normally distributed random numbers with zero mean and unit 
variance. Let the control action be turned on at the data point 18. In other words, we wait 
no time to close the system loop as soon as the first nonzero vector v p (k) defined in Eq. (21) 
is formed. Figure 1 shows the open-loop and closed-loop histories of input and output. The 
solid curve is the open-loop response and the dashed curve is the closed-loop response. The 
control gain starts with a zero vector and ends with the controller 

u c (k) = -0.3848 u c (k - 1) + 0.72llu c {k - 2) + 0.2536u c (A: - 3) 

— 0.0681u c (A: - 4) - 0.0150u c (A; - 5) + 0.0000u c (A; - 6) 

— 0.9828?/(A: - 1) - 0.7304y(fc - 2) - 1.2904y(fc - 3) 

-0.1818y(fc - 4) + 0.1908y(A; - 5) - 0.0007y(A; - 6) 

The controller converges to the one using the batch approach shown in Eq. (19) without any 

control action. 

Let the output be added with some measurement noise so that the signal to noise ratio 
is 11.4 (i.e., the output norm divided by the noise norm). The noise is random normally 
distributed with zero mean. Set the values of p and q to the same as that for deadbeat, i.e., 
p = q — 6. The open-loop and closed-loop time histories axe shown in Fig. 2. The control 
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Figure 2: Input and output time histories 
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gain starts with a zero vector and ends with the controller 

u c (k) = —0.1208 u c (k - 1) + 0.6969n c (A: - 2) - 0.02486u c (Jfc - 3) 

-0.0780u c (fc - 4) + 0.0704tt c (fc - 5) + 0.0762u c (fc - 6) 

— 0.5724y(A: - 1) - 0.3321 y{k - 2) - 0.7017 y(k - 3) 

+0.0593y(A: - 4) + 0.2075 y(k - 5) - 0.1076y(/c - 6) 

Although the hnal control gain looks considerably different from the one shown earlier for 

the noise-free case. The performance for both noise-free and noisy eases is quite similar. 

Increasing the value of q does not seem to improve the performance. In some eases with q 

larger than p , the performance is worse than the ease with p = q when the control action 

is turned on too early. Given sufficient time steps for the control gain to converge to a 

reasonable level, the performance may be improved somewhat particularly for the ease where 

q > p. One may raise the question whether the recursive controller design works for the ease 

where the order of the controller is smaller than the order of the system, i.e., p is smaller 

than 6 for this example. Let us choose p = q = 4. Using the same set of data for the previous 

ease, the input and output time histories arc shown in Fig. 3. Obviously, the performance 

is not as good as the ones shown earlier but it is still stable and somewhat acceptable. In 

practice, the order of a system is unknown and thus there is a great possibility that the 

values of p may be smaller than the true one. 

7 Concluding Remarks 

A new recursive predictive control method has been presented in this paper. System identi- 
fication was reformulated in such a way that it fits better for predictive controller designs. 
In other words, the conventional thinking in system identification has been re-oriented to 
focus on the control design process. The conventional procedure for any controller designs 
includes two steps, i.e., first perform system identification within an acceptable level of ac- 
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curacy and then conduct a controller design. The error in system identification will likely be 
accumulated and carried through the controller design process. As a result, the conventional 
approach tends to introduce more error in the controller design than system identification 
itself. The method derived in this paper uses a different approach for system identification 
to eliminate the additional controller design process. For noisc-frec cases, both conventional 
and new approaches produce an identical result if the predictive controller is unique. For 
non-unique controllers, the new approach provides the control gain smaller in norm than 
that from the conventional approach. This is due to the fact that, instead of minimizing 
the output error, the new approach minimizes the input error to compute the control gain. 
The proposed recursive method has a considerable advantage of computational speed. For 
noisy cases, numerical simulations have showed that the new method is more robust than 
the other methods. 
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